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ABSTRACT 


Relying heavily on the Elsasser and Culbertson (1960) computa- 
tional system, computations of terrestrial radiation incident at a 
black body interface have been programmed for an arbitrary atmosphe- 
ric sounding. The program has been applied to a random sampling of 
62 radiosondes from the "Wark-sounding-catalog." 

Since Wark et al. (1966) have kindly made available simulated 
Nimbus II channel 2 and channel 4 specific intensities for these same 
atmospheres, a multivariate regression was derived relating downward 
flux computations to Nimbus II channel 2 and channel 4 readouts as well 
as to two other gross air parameters: (a) total reduced water vapor 
depth, and (b) interface pressure. All four independent variables gave 
high statistical significance, with channel 4 filtered flux accounting 
for the major portion of the "explained variance" of the dependent 
variable. Total ozone was also tested but yielded no statistical sig- 


nificance on the regression. 
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1. Introduction 

A knowledge of downward infrared radiation flux is required in 
order to understand properly the heat partitioning in the earth-atmo- 
sphere and sun system. The atmosphere is much more transparent to 
short-wave than to long-wave radiation. This is depicted by the 
absorption spectra for the atmospheric gases on p. 153 of Fleagle and 
Businger [6]. From these spectra, one can see that there is strong 
absorption in the infrared by water vapor, carbon dioxide and to a 
lesser degree by ozone. It is the goal of this paper to develop a 
computer program that is readily adaptable to computation of downward 
radiation flux from soundings using input sounding data from each of 
these three atmospheric constituents. As a by-product, the secondary 
goal is to relate empirically the results to a regression equation of 
selected satellite air mass properties. 

Wark et al. [16] studied intensively the problem of upward flux 
and consequently compiled 106 randomly distributed radiosonde observa- 
tions from the I.G.Y. files in what is hereby termed the "Wark-catalog." 
This composite set of soundings cover a wide range of radiative condi- 
tions. Moreover, simulated Nimbus II data was kindly provided by J.H. 
Lienesch of the National Satellite Center for the same 106 soundings. 
Downward infrared flux measurements may be made by the Kuhn-Suomi 
radiometer but these give only limited areal coverage compared to the 
MRR data from satellites; thus, it was decided to compute downward 
flux in a manner similar (not identical) to that employed by Wark 
et al. in his 1962 study, and to perform a multiple linear regression 
analysis with certain simulated Nimbus II channel-readout values as 


input variables. The computation of downward infrared flux is made 
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applicable to the top of the undercast cloud or to the ground level, 
but the computer programming is left flexible so as to apply to down- 


ward flux at an arbitrary number of sounding levels above the interface. 


2. Nature of the Data 

The 106 radiosonde observations in the Wark-Gatalog are given as 
profile listings of pressure (mb), temperature (?K), mixing ratio of 
water vapor (g/kg), and mixing ratio of ozone (cm NTP per mb). Each 
sounding or case atmosphere is identified by a number, the name of the 
radiosonde station, and the date and time of the sounding. The sky 
conditions is indicated by "clear" or by the height (in mb) of the 
top of the undercast cloud. Temperatures and pressures are given at 
approximately 25 + 5 selected levels as well as at the interface. 
The interface is defined for this study as -€— the earth's surface 
for clear cases or the top surface of the undercast for cloudy cases. 
Total (corrected) water vapor and ozone masses are also given in this 
Wark- catalog. The models developed by the Meterorological Satellite 
Laboratory were taken from radiosonde data over all latitudes and 
seasons. All the soundings were from the I.G.Y record-file and went 
to at least 25 mb. Above this level, temperatures were extrapolated 
parallel to the appropriate Supplementary Standard atmosphere to the 
pressure of 0.1 mb. Furthermore, stratospheric humidity was extrapo- 
lated to conform in general to the Supplementary Standard distribution 
for the particular air mass at the station (standard cP, standard mT 
air masses, etc.). The 62 cases selected are composed of 47 undercast 
cases at various altitudes and 15 clear cases, the latter being ran- 


domly selected from the Wark-catalog. 
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Nimbus II satellite data is used to compute filtered channel 4 
upward radiation flux and the interface Black Body radiation flux for 
each of the 62 selected atmospheres. The Nimbus II data lists total 
and filtered channel 4 (5, to 304) infrared specific intensities for 
five zenith angles (theta), that is, б = 09, 20°, 45°, 60°, and 78.5°. 
The specific intensities for 9 = 90° are obtained by a Lagrangian 
form of extrapolation which is explained below in connection with 
equations (19) and (20). The radiation intensities of the 62 cases are 
listed in ergs/(cm2sec stdn) and are converted to watts/(m2 stdn) in 
this paper. It should be noted that the characteristics of the sat- 
ellite infrared sensing-filter system must be eliminated so that the 
results do not depend upon this optical system to any appreciable degree. 
The National Satellite Laboratory minimized the channel 2 and channel 
4 sensor-filter effects on the Nimbus II data by normalizing the filter- 
sensor response functions for the various IR wave channels, hence, 
the resulting response was the "best possible" estimate of the thermal 
flux (5-30) for the channel 4 radiometer, and of the black body flux 
(10-11u) or its equivalent black body temperature for the channel 2 
radiometer. Also, the sounding data lists temperature at the inter- 
face of each of the case atmospheres, which Keith [9] found to be very 


accurately described by the Nimbus II channel 2 readout. 


OY Computation of Downward Flux of Radiation 


Downward infrared flux of radiation arriving at the interface is 
quantitatively computed for each of the 62 selected atmospheres from 
the Wark-catalog by the Elsasser computational system. The downward 


flux computation is divided into five parts. 


T5 


E diaz ү pego HEUS S ОЧ (utm (1) 
The last two terms of equation (1) represent the "so called overlap 
of flux" and, hence, are subtracted from the first three terms of 
equation (1) which represents the downward infrared flux not cor- 
rected in any way from the manner described by Elsasser and Culbert- 
son [5] (henceforth denoted by EC for abreviation). The computation 
of the terms of equation (1) 15 explained below in connection with 
equations (6) through (11). 

In actuality, an energy-sharing process occurs in the 15y and 
9.6u bands which correspond to the CO2-H20 and O3-H20 overlap fluxes, 
respectively, thus, accounting for the AF-terms. For example, 
according to Hanel, Bandeen, and Conrath [7], ozone should be con- 
sidered the primary radiator in the 9.6y band, hence, water vapor 
is the weak secondary absorber-radiator resulting in the correction due 
to the shared radiation, AF(u*, yx). 

It is necessary to know the layer values of reduced optical path 
for water vapor, ozone, and carbon dioxide in the computation of radia- 
tive fluxes by each of the three constituents. The reduced optical 
path du* of an absorber has been taken from EC, and its relation to the 


true optical path is given by 


1 
dut = P (20) as 
FON T 


(0) 
According to EC [5], the dependence on temperature is slight and, hence, 


is excluded in this study. The finite difference relations for the 


resulting reduced optical paths of the three main constituents are given 


by 
K K 
ut = 2453 - E(qip (pj39 (&4445p). X 1.00705 X 1076 (g/cm2) (3) 
i=l i=l 
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К К 


0% - ZA,U* - (044) (44р) X 107 (cm NTP) (4) 
1=1 1=1 
K K 

U* = ДАМЕ - тө (А. ыр) X 2.4469 x 10-4 (cm NTP) (59 
1= = 


Basically, the EC computational technique for the determination of 
downward infrared radiation flux is realized by the integration of 
the area under the curve as shown in Fig. 1 with the use of Elsasser- 
Culbertson [5] flux tables of R(u*,T) which have been incorporated 

in the author's computer program RADFLUX found in Appendix I-A. In 
essence, one enters the EC tables with logarithmic values of the 
constituent optical paths and their corresponding temperatures, and 
obtains the respective "R" values which are explained below in connec- 
tion with equations (6) through (11). These "R" values are then 
integrated over the entire temperature range of the sounding result- 
ing in the flux of radiant energy due to the particular constituent. 
There are various forms that can be used to determine the energy 
integral but Elsasser et al. [5] and Yamamoto [17] used temperature 
instead of optical path (u) as the independent variable. 


The "R" term has been taken from EC and is given by 


R(u*,T) = 
Vo 


= (1-tp(u*)dy , (6) 
hence, one can readily see from equation (6) that infrared flux can 
be realized by an integration scheme of the R(u*,T) term over the 
temperature range. For example, downward infrared radiation flux 


for water vapor and its corresponding finite difference form are 


given by 
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Fig. 1. Schematic form of the distributi dat of R( UST) from tables in 
Elsasser and Culbertson 15) as applied to a typical sounding representative 
of the U-T distribution for any of the three atmospheric constituents. 

The area bounded by the sounding and the isopleths U* - 0 and U* = Uf 

is the downward flux uncorrected for possible overlap. 
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Те Т 


Fe 7 0) Rugs Dare  “°к[шҗ(т),т] ат (7) 
O Te 
і=іғ i=80°C 4 
Қт м tap Tia CIR ES E RigQ(u*g,T)) (T4,44-T) + 
= 1=1 
JE 


_ Buses. T)dr (8) 
For water vapor, the range of wave number (w,Vf) is listed, by 40 


cm i intervals, from 20 to 2600 спі, which is regarded as spanning 


the effective infrared range of the terrestrial radiation spectrum. 


The finite difference approximation for the Rj 4 1; in equation (8) 
assumes the form given by 
EX LL e(t ID c Rau (т, 15,13 24/1 /2 (9) 


and the last termof equation (8) is obtained from Table 20 of EC. A 
similar formulation for Foo, and Fo, is realized by the use of Tables 
11 and 13, respectively, in EC; however, Reg), and Коз values are re- 
stricted to the frequency ranges of 540-820 cm 1 and 970-1130 cm 1, 
respectively. 

The last two terms of equation (1), the H50-CO; and H20-03 over- 


lap flux terms, are defined by the following expressions 


T То — 
ДЕ (их ЛДК) = ( тр(ихр)Кс02 (ИК т)ат + $ tp(u*) Reo, (U*, T) dT (10) 
О T 
DN Б. i 
AF(u*,U*)= J tp(u*g) Ro4 (U*eg , T) dT + / Tp (wt) Rog (UF, T) aT (11) 
О | 


Е 


Moreover, values of Z (u*) in the H20-CO, overlap was determined by 


Martin [10] according to 


~ 


m Vf V£ 2 x= AL u* 2 
соу = (- С e-t dt) m (12) 
F FD ve T Я 
А № Vo үт, ` 


У 
> ВутАу vf 
Vo 5, BST Av 


Vo 
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where L is the generalized absorption coefficient listed by 40 ст“ 1 
intevals and Bp is the Planck specific black body flux for a given 
temperature. The results of this computation are shown in Table 1. 
The H20-CO2 overlap flux occurs in the 540-820 cm”? range, hence the 
water vapor TF(u) within this range was obtained by performing the 


summation indicated by equation (12) using both 1, , and BT at 20°C, 


FV 
and at 40 cm 1 intervals. L. values were taken from Table 10 in EC. 
According to EC, 20°C is a representative base temperature for т 
observations since the H50-CO y overlap is in practice most important 
in the lower layers of the atmosphere. In contrast, Wark et al. [16] 
used Yamamoto's [17] Труо Values which were based upon theoretical 
transition probabilities; however, Palmer's [14] laboratory measure- 
ments made at 300?K seem to give good agreement with Yamamoto's 


values for 25 cm 1 


intervals. After some additional smoothing, it 
was noted that the quoted L, values over 40 cml intervals in Table 
10 of EC corresponded closely to Yamamoto's quoted values of А2” 
реса ће of the slightly wider interval used in EC. 

Furthermore, 1F(u*) in the H20-03 overlap is given by 

Tp (u*) = exp[-0.1167 u*] (13) 
Equation (13) uses the continuum (970-1130 cm 1) absorption coef- 
ficient of water vapor proposed by Hanel et al. [7]. In both the 
H290-CO2 and H20-03 overlap regimes, the water-vapor sharing path has 
been increased by the factor 5/3 in order to transform beam trans- 
missivity to flux transmissivity. 

The results of the foregoing Elsasser computational scheme 


using the author's computer program RADFLUX on the 62 selected 


case atmospheres from the Wark-catalog are shown in Table 2 along 
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Table 1. 


Water vapor flux transmissivities in the H50-CO» overlap 


regime (540-820 cm-l). 


log u* те(ч%) 

-6.0 0.9999997150 
-5.7 0.9999994259 
09 0.999998527 
-5.0 0.999997150 
-4.7 0.999994259 
-4.3 0.99998527 
-4.0 0.99997150 
-3.7 0.99994259 
-3.3 0.9998527 
-3.0 0.9997150 


log u* 


50. З 


qs 7 
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тр(иж) 


0.9994259 
0.998527 
ОУ Шо О 
0.994259 
0.985665 
0.97128 
0.94431 
0.86680 
0.77328 
0.67567 


0.54833 


with their corresponding case-number identifiers from this catalog. 
The table lists uncorrected fluxes for water vapor, ozone, and car- 
bon dioxide together with the CO2 and ozone overlaps. The CO» and 

ozone overlaps are equally partitioned arbitrarily between CO, and 

water vapor, and ozone and water vapor, respectively. As a result, 
half of the overlapped flux is restored to"Cosamd ezone, thak is, 

the partitioned overlap fluxes are subtracted from their respective 
uncorrected constituent fluxes as shown below in connection with 


equations (14) through (16), thus, realizing corrected constituent 


fluxes. 
Fwy (CORR) = Fwy - AF(u*,U*)/2 - AF (u*,4*)/2 (14) 
Fco2 (CORR) = Есо2 - АЕ(а* М*) /2 (15) 
Fo3(CORR) = Fo3 - AF(u*,U*)/2 в (16) 


It should be noted that all the fluxes shown in Table 2 are in watts/ 
m? and that the mean of each overlap and constituent flux as well 


as the average percent contribution of the corrected constituent 


flux on the total downward flux are given at the bottom of the table. 
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Table 2. Selected 62 case atmospheres as indentified from the Wark-catalog, and the computed 


downward flux from program RADFLUX. 
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4. Determination of Channel 4 Upward Flux of Radiation. 


For each sounding case, the filter-normalized upward flux in the 
wavelength of 5u to 30u was available as a result of computations fur- 
nished by Lienesch. This was in the form of the "so-called" channel 
4 specific intensity at the top of the atmosphere. Channel 4 specific 
intensity I, (99 was then transformed to a channel 4 flux for each of 
the selected 62 cases by using the trapezoidal rule for integration on 
I4(8) as given by 


ТГ 
$F, = 2m /2 I4(0) d(sin^o) (17) 
О 2 


and its corresponding finite difference form is given by 


/2 
oF, = 21 2 I, (8, )+1, (6441) 


0 fim 7 Х(ѕіп2(0;41} - sin?(04)) (18) 


where í = 1, 2, 3..5, and 01 = 09, 0> = 200, Өз = 45°, Ө„ = 60°, 

5 = 9-3, апа P 909. The trapezoidal rule was selected instead 

of other finite integration schemes because of the unequal zenith- 

angle intervals of integration. Also the channel 4 upward flux relation 
is the same as that used by Wark et al. [16] which allows for the usual 
approximation of azimuthal symmetry. 

The zenith angles (6) of the Nimbus II listings went as far as 
78.5°, hence, the specific intensity for @= 90° was extrapolated by a 
Lagrangian form expressed by 

I(0) = I1(0-02) (6-03) + I7(6-6,) (6-63) + 13(0-01) (6-69) (19) 

opto) (82-8) (862-83) — (53-61) (63-69) 
and which reduces to 


1(900) = 178. 59 + 0.68661450 = 1.86491600 "P 1.1783178 590 (20) 
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According to Wark et al. [16], radiation for angles greater than 
78.59 contributes only about 4 percent of the total flux and the energy 
obtained by integrating eq. (17) between O = 78.59 and 0 - 90° lies 
well within this 4 percent range. 

The black body radiation flux (от) at the interface is given 


directly by the interface temperature (°K), and when o= 0.56687 Х 10-7 


4 


В 15 1п watts/m?. 


(watts/m- ок“), oT 
The two satellite air mass properties, oF, and om. are made avail- 

able for this study by the foregoing procedures, where T, is taken 

equivalent to the channel 2 emittance temperature in a simulated 

Nimbus readout. Moreover, the foregoing integration employing eqs. 

(18) and (20) were incorporated in the author's computer program SUMFLUX 


(Appendix I C). The 62 selected Nimbus II cases as well as the computed 


air mass properties are shown in Table 3. 
5. Statistical Results and Inferences 


The object in the statistical treatment which follows is to obtain 
a specification equation for Fg in terms of satellite readouts and 
other readily accessible air-mass data. To accomplish this, there was 
available a BMDO3R computer program that executes the multiple regres- 
sion analysis establishing a regression plane given by the linear form 
of 

| до + A124 + А2) + А32. % А,2, % Ас26 (21) 
with the final results of the regression plane shown in Table 4. 

The original set of 62 independent data atmospheres was split 
randomly into two equal sets. One set, Sample A which is denoted by 


the asterisks in Table 3, was used to develop the regression while the 
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Channel 4 filtered flux and interface black body flux on 


trable 2. 


the 62 case simulated Nimbus II data. 
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remaining set, Sample B, was used to test the stability of the test 
data on the check-data which was also derived by the author's computer 
program RADFLUX. 

According to Miller 13] and Martin et al. [н], the regression 
screening analysis uses a F' critical value both to accept and remove 
individual variables from the regression. The critical F' value uti- 


"probability level" that is a function of the number "k" of 


lizes a 
predictors tested and of the Sth selected prediction system. The pro- 
posed critical value corresponds to the tabulated F' value at a prob- 
ability level of a=a*/(k-S+1) with 1 and N-S-1 degrees of freedom, 
where N is the sample size to be tested. Since the predictors (k=5) 
have been pre-selected, the selection choice "S" is equal to 1, hence 
K-S+1 is equal to 5. 

The predictors employed are ¢F4, ae už, the interface pressure 
"р", апа U* which correspond to Z1, 22, 23, Zą and Zs, respectively, in 
eq. (21). The order of the predictors is a vital consideration in 
describing the "explained sum of squares by the regression." Channel 
4 filtered flux ($F4) is referred to as the thermal channel, for its 
infrared range is 5y- 30y, whereas channel 2 filtered flux, the window 
channel, has an infrared range о 10-111. Channel 2 filtered flux, 
however, is equivalent to the interface black body flux, cT£ ,and there- 
fore contains very little atmospheric structural information.  Pre- 
scribing $F, ahead of от allows $F, to assume a larger proportion of 
the explained variance than oT}. 

The most readily adaptable test for the significance of the 


predictor is the F test based upon an analysis of variance. The F 


statistic is defined as 
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Е = Mean squares explained by the predictor(s) (22) 


Mean squares not explained by the predictor(s) 

The F statistics derived from the results of the BMDO3R program are 
shown in Tables 4 and 5. These statistics are used to determine the 
significant specification on F4 by the multivariate regression and by 
each of the predictors. Significance in this analysis means that the 
independent variable and/or variables selected for the F test did not 
affect the outcome of the regression by mere chance at the specified 
level of belief. For this analysis, an "a" of 0.01 is chosen resulting 
in a composite 99% confidence level of belief for the multiple regression. 

Also shown in Tables 4 and 5 are the tabulated critical F' values 
and the "o" critical level for rejection of a possible air-mass variable 
together with the corresponding "F-value upon entry." Since "ax"! 
is chosen to be 0.01, the "o" critical levels for the 5 predictor 


regression are 0.002. The "F-values upon entry" are defined by 


F(1,N-m-1) = Mean Squares explained by the regression (23) 


Mean Squares not explained by the regression 
where "N" denotes the number of cases in the sample, "m" is the 
sequential order of the predictors (m = 1, 2....k), "k" is the number 
of predictors in the regression, and the number pair (1,N-m-1) denotes 
the number of degrees of freedom associated with the numerator and the 


denominator, respectively. For example, the 'F-value upon entry" for 


the first predictor ($F,) in the regression is given by 


(24) 
F(1,29) = (Sum of Squares attributable to the regression)/l 


(Sum of Squares not attributable to the regression) /29 
The addition of the next predictor ,oT¢, requires the definition of 
a new F-statistic specifically restricted to the testing of ote, 
after the contribution of $F, has been excluded. This new F-statistic 


follows the same form as given by eqs. (22) and (23) but with one less 
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degree of freedom in the denominator. The remaining independent 
variables are examined successively in the same manner with a new 
F-statistic specifically restricted to the testing of the selected 
variable after the contribution of the variables already tested have 
been excluded. 

The composite F-statistic for the multiple regression is based 


upon 


F(k,N-k-1) = (Sum of Squares attributable to the regression) /k (25) 


(Sum of Squares not attributable to the regression) /(N-k-1) 
Use of eq. (25) with the results from Table 4 gives the F-statistic 
for the 5 predictor multiple regression as F = 116.03 which is far in 
excess of a critical F'(5,25) = 3.86 appropriate at a 99% confidence 
level. 

The results of Table 4 show not only that the multivariate 
regression, but also that each of the single variate predictors are 
significant with the exception of the last independent variable, U*. 

The "F-statistic upon entry" for y* is 0.90416 in contrast to the 
required (tabulated) critical F'-value of 14.0 at 99.8% confidence level 
of belief. This indicates that U* could only have affected the outcome 
of the multivariate regression by chance at a very low level of probab- 
ility. An explanation for the lack of significance by ozone could lie 

in that the computed amounts of ozone in a column of atmosphere are near- 
ly constant across the 62 case sample resulting in a negligible influence 
on the regression. Since U* showed so little significance in the 
regression analysis, another analysis was made with U* excluded. The 
regression plane for the 4 predictor case is given by the linear 
expression 


Y' = Ag + AiZ] +A 


1 4 аа AN (26) 


2 379 4 4 
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with the final results of the multivariate regression shown in 
Table 5. The same procedures and testing criteria as before were 
used resulting in each of the 4 predictors as well as the composite 
set of predictors showing significant contribution on the regression, 
thus, the 4 predictor regression was selected as the testing regres- 
sion plane for the independent data sample. The five variate regres- 
sion equation is given by 

F4 = -85.535 + 1.750(9F,¿) - 0.486 (от8) + 47.466 (и*) + 0.129(р) + 

-121.866(U*) (27) 

and the four variate regression equation is given by 

Ба = -128.089 + 2.022(9F4) - 0.605(078) + 49.897(u*) + 0.129(p) (28) 
It is interesting to note that the multiple correlation coefficients 
shown in Tables 4 and 5 decrease very slightly from 0.9791 in the five 
variate case to 0.9784 in the four variate case. An explanation for 
this could lie according to Panofsky et al. 5] and Martin et al. [11] 
that the longest equation may have actually overfitted the first analysis 
ascribing some of the variation due to small-scale perturbation to U* 
by chance. 

Finally, the 4 predictor regression was tested using Sample B 
data by an analysis of variance or of mean squares, that is, a test 
of variation about the regression plane. Again, significance of the 
regression plane was tested by the F test as given by equation (24). 


Total sum of squares with its components is given by 


N LA N р М a 
r бук =: Cy uoc Er у (у.-у) (29) 
i=l ^ PE i-l ^ 

(A) (B) (C) 


where VT denotes the regression estimates for F4,"y," denotes the 
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observed F4 values as determined by the computer program RADFLUX, and "N" 
denotes the total number of cases in the independent data sample. Term 
(A) gives the total sum of squares; term (B) gives the sum of squares 
of deviations from the regression plane or the sum of squares not 
explained by the regression; and, term (C) gives the sum of squares of 
deviations of the regression estimates from the mean (у) of the observed 
values or the sum of squares explained by the regression. 

The statistics on the independent data sample are given in Table 
6. The required (tabulated) critical F'value [Р' (4,26)] is equal to 4.14 
at 99% confidence level of belief. Опе can readily see that the "Fr" 
statistic for the independent data F(4,26) = 118.2546, is considerably 
larger than its corresponding tabulated critical F'value, thus the 4 
predictor regression verifies as being significant for determining 
F¿ within 99% confidence level of belief. 

Furthermore, the results of the 4 predictor regression on the 
independent data sample were tested for shrinkage using the relation 
given by 


~ 


Ry? - l - (Sum of Squares not explained by the regression) (30) 
(Total Sum of Squares) 


where Ry is the effective coefficient of determination and Rk is termed 
the "effective multiple correlation coefficient." The results are shown 
in Table 6. 


The percent shrinkage between the explained sum of squares of 
samples A and B is given by RA - A. Here й is the coefficient of 
determination for eq. (28) (Sample A data). Table 6 shows that the 
regression performed exceptionally well on the independent data sample 


2 


with only a 5.6275 amount of variance shrinkage and with a R4 and Ry 


of 90.1% and 95.1%, respectively. 
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Table 6. Analysis of variance on the independent data (Sample B) 


Degrees Sum 
Source of Variation of of Variance 
Freedom Squares 
Explained Variance 4 86065.300 2868 .8430 
Unexplained Variance 26 9462.072 57 6802 
Total Variance 30 955277372 302075327 
F-value Critical Coeff. of Multiple Percent 
F' value Determina- Correl. 


F(4,26)  F'(4,26) tion (R42)  Coeff. (R4) Shrinkage 


118.2546 4.14 0.901 0.951 5.02 
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The foregoing analysis showed that the 4 predictor regression was 
indeed verified on the independent data sample; but, the process 
of randomly dividing the original 62 Nimbus II cases in order to realize 
this verification excluded 31 perfectly good cases from the regression 
analysis, thus, it was decided to obtain a regression based on the 
entire 62 Nimbus II cases. One should anticipate a more representative 
regression than the one previously derived, because the sample size is 
twice as large as that used in the original analysis. Although 
the new analysis used the same methods and testing criteria as that 
used in the 31 case sample, it was, however, confined to the 4 predictor 
analysis because of the earlier findings in the study. The final regres- 
sion is given by 

F} = -129.847 4 1.824(4F4) - 0.391(oT2) + 43.293(u*) + 0.097(р) (31) 
and its corresponding statistics are given in Table 7. It is interest- 
ing to note that the multiple correlation coefficient decreased slightly 
in going from the 31 case analysis to the 62 case analysis. An explan- 
ation for this must lie in that the sample B "explained variance" was 


slightly smaller than that for sample A. 


6. Conclusions 

The programmed computational system for downward infrared radiation 
flux due to the absorbing/emitting properties of water vapor, ozone and 
carbon dioxide is given in Appendix I A by the author's computer program 
RADFLUX. This program may be readily used for the determination of Fa. 
The only necessary inputs are pressure (mb), temperature (°K), water 
vapor mixing ratio (g/kg), and ozone mixing ratio (cm NTP per mb) for 


selected levels in a sounding. Moreover, this program can' be easily 
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modified as shown in Appendix I B in order to determine Fa from any pre- 
selected level to the top of the atmosphere (0.1 mb). It should be 
noted that the mean ratio of downward radiation flux due to water vapor 
(Ps? with and without overlap corrections to the total downward 
radiation flux (Ра) was /3.9/7% and 85.6%, respectively. 

Furthermore, total and filtered channel 4 upward radiation flux 
can be determined by using the author's computer program SUMFLUX found 
in Appendix I C. Also, the interface black body flux or its channel 2 
equivalent temperature can be readily determined using this same pro- 
gram. 

Finally, the multivariate regression equation developed using 
Sample A Nimbus II air-mass properties for the determination of F. 
had an accuracy of 95.1% on the independent test data, Sample B. This 
can be seen in the foregoing analyses. The explained standard devia- 
tion of the Sample A regression was 95.84%, but the resulting explained 
standard deviation on the independent sample was 95.1%, giving a loss 
of accuracy of 2.74%. A partial explanation for the loss of accuracy 
lies in the process itself of developing a regression on a finite ran- 
dom sample. One should not anticipate a regression determined from a 
finite sample to perform exactly the same way on another sample of the 
same population. This is inherent in the randomness of the selection 
process and of the variations within the population. This effect can 
be anticipated to be partly minimized by using the largest possible 
sample. Another partial explanation could lie in the assumption that 


the regression is linear, but in actuality, F, might vary non-linearly 


d 
with one or any combination of its predictors. This would indicate 


that the true regression is some sort of surface and not a plane, 
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however, since the shrinkage of the multiple correlation was only 
2.74%, the regression plane appears to be a good estimate of the true 
surface or of physical reality. With the advent of improved satellite 
sensing systems, the regression equations can be expected to be 
improved by an analysis similar to the one given in this paper. 
Furthermore, an improvement on the regression plane may be anticipated, 
if the satellite data as well as the radiosonde data could be strati- 
fied to a common level, namely that of the interface. This, of course, 
would exclude pressure as a possible air-mass predictor, however, one 
could reasonably expect this stratification to reduce the scatter about 
the regression plane. Moreover, since the effect of ozone in a column 


of atmosphere gives a very nearly constant contribution to F, and, 


d 
thus, has very little effect on the regression plane, and since one 
should expect that the regression based on the largest possible sample 
to be the most representative to the true regression surface, the 


best regression for the determination of F, is, therefore, the 4 


d 


predictor regression based upon the 62 case Nimbus sample, eq.: (31). 
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Appendix I, Section A. Computer Program RADFLUX (FORTRAN 63) 


-СООР»ВОХ Т» ТОРА2 J»S/15/25»609400000 
-ЕТМ9(»Е» 


С 


60 


61 


59 


101 


103 


PROGRAM RADFLUX 

PURPOSE OF THIS STUDY IS TO EXAMINE DOWNWARD FLUX 
ODIMENSION TT(13)+.U(21) ».VUU(16) SUUU(L 23) sWATER( 21913) sOZONE( 16913)» 
1002 (23913) sRWATER( 21) »sROZ(16) »RCEO2(23) »P (40) sPMEAN( 40) »PDIFF(^0)» 
2T(4 )9Q(40) sQMEAN( 40) sQDIFF (40) »CQ( 40) sCQMEAN(40) sCQDIFF(40) » 
3SFW(40) »SFO(40) »SFC( 40) »CFW(40) »CFO( 40) » CFC( 40) sFUT( 40) »GUT (40) » 
4DUT(40) sFF(1)2GG(1)s0D(1) sCOTAUF( 21) sO3 TAUF( 40) sRCOLAP(1) > 

5СО21 АР (21) sO3RLAP( 40) sRO3BLAP (1) 2UC( 21) »sCOXMIS(21) sTCOLAP(1) > 
6TO3LAP(1) 

DIMENSION KASE(40) »sLEVEL( 40) 

DIMENSION Z2(40) 9S(40) »V(40) oR( 40) 
ODIMENSION RH20(40) »RO3(40) ,КВСО (40 ) ›ТЕМР]1 (40 ) , ТЕМР2 (40) »ТЕМРЗ(40), 
1 КОРЕ (13) я 5ОРЕ( 13) » ТОРЕ { 13 ) эТТМ ( 13) о 1710( 13) ›ТТС( 13) »›РМ( 13) › РО ( 13) » 
2RC113) »RURF (13) » SURF( 132» TURF (C13) 

DIMENSION MKASE(100) »ZCOLAP (100) »ZO3LAP (100) »2 FWV(100)»WFWV( 100)» 
14F03(100)»wFO3(100)»ZFCO2(100) ,WFCO2(100) »ZTOTAL(100) 

DIMENSION RATIO1(100) »sRATIO2Z(100) sRATIO3(100) 

КЕА0(50+60) ((МАТЕК(1»Ј)•» Ј=19»]13)»1=1•»21) 

КЕАр(50»60) ((С02 (1+ 0) + Ј=1»+13) ,]=1+23) 

КЕА0р(50»60) ((О2ОНЕ (19) ) • )Ј=2 1913) +1=1» 16) 

FORMAT ( 13F643) 

READ(50>61)(UlI)»I=1>21) 

КЕАО(50»61) (М( 1І)»1:1»23) 

READ(50»61) (UU(I)»1721»16) 

РОКМАТ (12284 6 О» ЗЕ36 0, 22 60+,723%0) 

READ(50»59) (RWATER(K) »K21521) 

READ(50»59) (RCO2(K)5»Kz1»223) 

READ(50959) (ROZ(K) sK=1916) 

FORMAT (12F6e2) 

READ (505701) (COXMIS(I)»I21521) 

FORMAT(5F15410) 

READ(50»7023) (UC(1)»121»521) 

FORMAT(16F4.2) 

DO 62 11913 
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62 
100 


63 


64 


777 
1000 


7000 
10 


ТТ(І)-1%10-90 

WRITE( 515100) 

FORMAT ( 1H1) 

WRITE(51>63)(TT(1)>»1=1,13) 

FORMAT (4X»13F8.0) 

WRITE(51»100) 

WRITE(51»64) (ULII)»(WATERL1,J),J=1»13),1=1»2]) 
FORMAT(FS.1»13F8.3) 

WRITE(51»100) 
WRITE(51964) (UU(T) » (OZONE (I 9J) 5 J=1913) s1=1916) 
WRITE( 519100) 
WRITE(51264) (UUUCT) +" (СО2 (1, )) , )=21» 13) +» 1=1%»23) 
INTER=63 

KUT = 0 

EONST=0 24838 

READ (5051000) KZ 

РОВМАТ ( 14) 

READ (5057000) KOUNT 

FORMAT( 14) 

READ (5091) (KASECI) sLEVEL(CI) PCI) » Т( 1)» О(1) » СО (1) +» 1=1%К2) 
FORMAT(213,F9.1,F5.0,F12.4,F10.3) 

WRITE (51»100) 

WRITE (51»2) (KASE(I)>LEVEL(1),P(1)>»Q(1)>»C0(I)>I=1>KZ2Z) 


20FORMAT (1X>4HCASE +, 4X»5HLEVEL » 4x ¿BHPRESSURE » 6X » 3HH20, 7X3 5HOZONE//(1X 


1914»918»ЭҒ14.4Ғ12.84»С51123)) 


KT= KZ - 1 

SUMW=0 e 

SUMO=00 

SUMCzO, 
CONST1=U.4342944819 

DO ¿€ I=-1s»KT 

РМЕАМ( 1) = (Р(1)+Р (1+1))/2• 
QMEAN( I) =(0О(1)+0(1+1) )/2• 
СОМЕАМ(1) = (CQ(I)*CQUI*1))/2, 
РОТЕР(1) = Р(1)-Р(1+1) 
&DIFF(I) = Q(1)-Q¢C141) 
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28 
26 


54 
22 


48 
46 


21 


53 


47 
20 


» 


CQDIFF(I) = CQ(I)-CQ(I+1) 
THE PURPOSE OF THES SUB-PROGRAM IS TO CALCULATE THE LOG VALUES FOR 
THE REDUCED OPTICAL PATHS. 
SUMW=SUMW+QMEAN (1 ) ХРМЕАМ( 1) ХРРЈЕРЕ (1) #( 16 00705Е-6) 
SUMO=SUMO+CQMEAN(I)*PDIFF(I1)*(].E-5) 
SUMC=SUMC+PMEAN(I)%*PDIFF(I)*(2.4469E-4) 


SFW(1)=SUMW 
SFOLI)=SUMO 
SFC(I)2SUMC 


IF (SFW(I)) 56957558 
CFW(1)=CONSTIX*LOGFISFW(I)) 


CONTINUE 


IF (SFO(I)) 529153994 
CFOLI)=CONSTI*LOGF (SFOLI)) 


CONTINUE 


IF (SFC(1)) 46947948 
CFC(1) =CONST 1*LOGF (SFC(I1)) 


CONTINUE 
GO TO 20 
CFW(1)=19 
GO TO 56 
EFO(1)=19 
GO TO 52 
CFCt1)=19. 
CONTINUE 


OWRITE(51s4) (15Т(Г)› РЕГ) ›ОЯ (ТГ) › РМЕАМ (1) ›РОГЕЕ (1) › ОМЕАМ (|) › $ЕМ (|) »› 
ЈУРО(1)»БЕС( 1)» 1=]•К1) 
4OFORMATL2H 1I»4X>»1HT>6X91HP»>7X>1H0»4X»S5HPMEAN» 3X>»5HPDIFF+3X» 
15HQMEAN»6X93HSFW»6X» 3HSFO96X>3HSFC//(13>»F6.0»FB8.1»F7.3»2F8.1»4F9.4 


2)) 
WRITE( 519100) 


WRITE( 5109) (І9Т( +1). СЕМ(І) sCFOC(I) »CFC( 1) sIT=19KT) 
OFORMAT (3H Іэ2Х»ӨНТ(ІЗ1)»Э7Х»ЭЭНСҒЫ»17Х»ЗНСҒО»917Х»9ЭНСҒС//( 15»Ғ5.0» 


LER 5) 


NOTEee eses ANY LOG VALUE DENOTED AS 19 INDICATES THAT IT HAS BEEN 


FORCED OUT CF 


RANGE OF THE 


R 


VALUE TABLES BECAUSE OF THE 


ac 


U 


VA- 


OO 


89 


LUE BEING 2ЕКО. А(50» NOTE THAT U VALUES CAN NOT BE NEGATIVE 
DUE TO THE ARITHMETIC AND SEQUENTIAL PROCESSES OF THE PROGRAM. 
WRITE( 512100) 

JUMP= 1 

KH= 

KO= 

KC= 

THE PURPOSE OF THIS SUB-PROGRAM IS TO DETERMINE THE SPECIFIC 
WATER R VALUES FOR EACH LEVEL OF EVERY CASE, 

DO 13 I=KOUNT sKT 

TÜ = T(I+1) — 273. 

WU=CFW(1) 

DO 23 K=1l9+12 

IF(TU-TT(K)) 232940950 

IF(TU-TT(K+1)) 40943423 

CONTINUE 

IF (JUMP )289 9189989 

WRITE(51288) KASE(C |1) Т( ІЗІ) әм) 


880FORMAT(1X»8H CASE 7 I4»4X»8H TUH2O -Ғ5.0»5Х911Н LOG UH20 sF6435»23 


133 
40 


43 


191 


80 
39 


10 
74 
73 


lH ARE OUT OF TABLE RANGE) 


IF (JUMP 10726» 13 

I T2K 

GO TO 65 

IT=K+1 

GO TO 65 

DO 39 L=1»20 
IF(WU-U(L))39+709+80 
IF (WU-U(L41)) 70273539 
CONTINUE 

GO TO 89 

Ім = L 

IF (JUMP)120»110»90 
Iw = L+1 | 
GO TO 74 


9O0ODCALL FFIN(WUsU( IW)»9U( IW+1)9»TUITTC(IT)9TTCIT+1)9 WATER IWo IT)» 


IWATER(IWs IT+L I)» WATER IW+ 1» IT) ,WATERLIW+1L»IT+ 1)» FUTLI)) 
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KH=KH+1 
RH2O(KH) = FUT(1) 
TEMP1(KH)=T(1+1) 
13 CONTINUE 
WwRITE(515»100) 
JUMP = JUMP - 1 
DO 26 I =KOUNT oKT 
ТІӘ-Т(1»1)-273. 
WU=CFO(1) 
бо ТО 36 
189 WRITE( 519188) KASE(I) »s TC I+1) swu 
1880FORMAT(1X»s8H CASE г 14» 4Хэ8Н Т0ОЭ = F5e095X910H LOG VOR =F6.39+23H 
1JARE OUT OF TABLE RANGE) 
СО ТО 133 
102 BO 105 L=:1,15 
IF (WU=UU(L))105»70»106 
106 IF(WU-UU(L*1)0)0705»73»105 
105 CONTINUE 
GO TO 189 
llOULCALL FFOZIUWUSUUCIW)»UU(CIW*l)sTUSTT(OOITOs TTCIT*1)»s0OZONE(C IW» IT » 
109Z2O0NECIWs.SIT*1)»sOZONE(CIW*1»IT) ;jOZONE(CIW*1» IT*1) GUT(CI)) 
KO=KO+1 
RO3(KO)=GUT(I) 
ТЕМРг(КО) = Т(1+1) 
26 CONTINUE 
JUMP = JUMP- 1 
DO 107 I=KOUNT KT 
TU=T(1+1)-2730 
WU=CFC(I) 
60 TO 36 
289 я«КІТЕ(51»288) КАСЕ(1) ЭТ(1%1)»эуд) 
288 FORMAT(1X»8H CASE = I4»s4X»8H TUCO2 = F54,0»5X»11H LOG UCO2 =F6.39+23 
1H ARE OUT OF TABLE RANGE) 
GO TO. 133 
103 bO 108 і:1,22 
IF (WwU-UUU(L)) 1085705109 
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109 IF (WU-UUU(L+1)) 702735108 
108 CONTINUE 
GO TO 289 
120UCALL FFCO2 (WU + UUU ( IW) sUUU( IW+1) sTUsTT( IT) os TTCOIT*10»5CO2€0 IW» IT)» 
LCOZCIWsIT+1) sCO2( I W+lsIT) sCO2( [w+elsIT+1) sO0UT(I1)) 
KC=KC+1 
RRCO(KC)=DUT (1) 
TEMP3(KC) = T(I+1) 
107 CONTINUE 
60 TO 77 
65 ТЕ(ЈОМР )103»+102+ 101 
THIS SUB-PROGRAM DETERMINES THE FLUX INTEGRAL FOR WATER VAPOR 
BETWEEN -80 TO -273 DEGREES CENTRIGADE. 
77 WU * CFW(KT) 
IND = 3 
00 150 Іг1»20 
IF (WU - U(ID) 15051705180 
180 IF(WU-ULI*1)0) 170»1735150 
150 CONTINUE 
170 Імгі 
174 GO TO (22092109190) IND 
173 IW=I+1 
GO TO 174 
190 CALL RRWAT(WUsU( IW) UC IW+l) sRWATER( IW) sRWATER( I W+1) s FF) 
WRITE (519191) FF 
191 FORMAT(1X»6H FF = E1729) 
THIS SUBPROGRAM DETERMINES THE FLUX INTEGRAL FOR OZONE BETWEEN 
-80 TO -273 DEGREES CENTRIGADE, 
INO*INO-1 
WU=CFO(KT) 
ВО 160 1=1»15 
IF (WU-UU((I))16051705181 
181 IF(WU-UU(LI*1))01705»1735160 
160 CONTINUE 
216 CALL RROZ(WUsSUU(IW)9*UUCIW*1)*ROZUIW) 9sROZUIW*1)9»GG) 
WRITE (519211) GG 
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182 
200 
220 


221 


400 


413 


40] 


414 


402 


FORMAT(1X»6H GG = El7.9) 

THIS SUB-PROGRAM DETERMINES THE FLUX INTEGRAL FOR CO2 BETWEEN 
-80 TO -273 DEGREES CENTRIGADE, 
IND=IND-1 

WUSCFCIKT) 

DO 200 Iz1»22 

IF (WU-UUU(1)) 20021705182 

IF (WU-UUU(1+1)) 17051735200 

CONTINUE 

CALL RRCO2(WUsUUU(IW) sUUUCIW*1) »*RCO2(IW) sRCO2(1W*1)0 500) 
WRITE (515221) DD 

FORMAT(1Xs6H DD = E1729) 

THIS SUB-PROGRAM INTEGRATES THE INTEGRAL ROT BETWEEN T INITIAL 
TO T FINAL IN DEGREES СЕМТКТСАОЕ о 
IND=IND-1 

АВЕА1=0, 

АКЕА2=0 

AREA3=00 

ТОТА(=0 • 

KHHzKH-1 

DO 400 I31,KHH 

У=ТЕМР1 (1) - ТЕМР1 (1+1) 
AREA]z=AREA1+(RH201 1)+RH2011+1))/20*Y 
WRITE (51»413) AREAI1 

FORMAT(1X»10H AREAL = Е17,9) 
KOOsKO-1 

BO 401 I215»KOO 

У=ТЕМР2 (1) - ТЕМР2 (1+1) 
АКЕА2=АКЕА2 + ( КОЗ(Т)+КОЗ (1+1) ) /2.*У 
WRITE (519414) AREA2 

FORMAT(1X»10H AREA2 = E1729) 
КССгКС-1 

ро 402 Is1»KCC 

YzTEMP3(I) - TEMP3(I-*1) 
АКЕАЗ=АКЕАЗ+ ( ККСО (1) +ККСО ( 1+1) ) Уго *У 
WRITE (519415) AREA3 
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С 
С 


415 


909 
90:3 


906 
914 
913 


981 


900 


924 
918 


982 


915 


FORMAT(1X»10H AREA3 = E1729) 
THIS SUB-PROGRAM INTEGRATES THE INTEGRAL ROT 
ТО 80 DEGREES CENTRIGADE. 
МС = 

МО = 

МИ = 

KUMP 3 1 

WUZCFWIKT) 

DO 900 L=1»13 

DO 903 J=1»20 

IF (WU-U(J)) 90329063909 

IF (WwU-U(J+1))9069973 9903 

CONT INUE 

60 TO 900 

Iw=J 

IF(KUMP) 983,9825,981 

IW=J+1 

GO TO 974 


BETWEEN 


Т 


F INAL 


CALL FFFN(WUSULIW)SsUCIW*1) »5WATERCIWSsL) sWATERCIW*19»L) »ROPE(L)) 


MW=MW+1 

TTW(MW) = TT(L) 

Rw(Mw) = ROPE (L) 
CONTINUE 

KUMP =KUMP—1 

WU=CFO(KT) 

DO 915 L=1 13 

BO 918 Jz1»15 

IF (WU-UU(J)) 91899069924 
IF (WU-UU( J4+1)) 90699735918 
CONTINUE 

Go TO 915 


CALL FFFNUWUSUUCLIW)SsUULIW*1) 95020NECIW9* L) 95020ONECIW*1»L) s SOPE(GL ) ) 


МО=МО+ 1 
TTO(MO)*TT(L) 
RO(MO) = SOPE(L) 
CONTINUE 


49 


КОМР =КОМР- 1 
WU=CFC(KT) 
РО 930 (1•13 
00 933 21,22 
IF (wWU-UUU(J)) 93359069936 
936 IF(WU-ULU(J+1)) 90699739933 
933 CONTINUE 
GO TO 930 
983 CALL FFFN(WUSUUU(CIW) sUUUCIW*1)  CO2( IW*3 LL) » CO2C IW 19 С) 3 TOPECL ) ) 
MC=MC+1 
TTC(MC) = TT(L) 
RC(MC) = TOPE(L) 
930 CONTINUE 
МОМР =] 
ТО=ТЕМР] (КН) -273. 
KL=KH 
KM=KO 
КМЕКС 
830 DO 800 K=1 12 
ІЕ(ТО-ТТ(К)) 800»803,806 
806 IF(TU-TT(K*1)0) 803,873,800 
800 CONTINUE 
IF (MUMP) 81558125809 
809 wRITE(512818) TU 
818 FORMAT(10H TF H20 * F54,0»3X»22H IS OUT OF TABLE RANGE//) 
KL=KL-1 
TU=TEMP1(KL)-2736 
GO TO 830 
803 IT=K 
874 IF(MUMP) 8835882881 
873 IT=K+l 
GO TO 874 
881 CALL FORGE](TUsTTIW( IT) sRW( IT) 4 RH20(KL) sAREA11) 
ADD1=0. 
IM=IT-1 
DO 824 L=1»IM 
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nn 


CALL FORGE2(RW(L)»RwW(L+1)»TTW(L)»TTW(L+1)»RURFAIL)) 
824 CONTINUE 
DO 827 L=1l>»IM 
827 ADD1=ADDI+RURFIL) 
THIS SUB-PROGRAM DETERMINES THE TOTAL FLUX FOR WATER VAPOR WITH- 
GUT CO2 AND OZONE OVERLAP CORRECTIONS. 
845 АКЕА1 = AREA] + ADD1 + AREA11 
WRITE (519410) ADD1sAREA11 sAREAI] 
410 FORMAT(1X»8H ADD1 ғ Е17%9»95Х210Н АКЕА11 * E174955X9»14HTOTAL AREAL 
1= Е17.9) 
MUMP=MUMP-1 
TU=TEMP2(KO) -273e 
GO TO 830 
812 wRITE(515833) TU 
833 FORMAT(10H TF 03 5 F5,0»3X»22H IS OUT OF TABLE RANGE//) 
KM=KM- 1 
ТО-ТЕМР2(ҚМ) -273. 
GO TO 830 
882 CALL FORGELI(TUsTTO(tIT) sROCIT) ;ROS3(KM)sAREA22) 
ADD2=0. 
IM=IT-1 
DO 836 [=1»1М 
CALL FORGE2(RO(L)»ROCL*1)»TTO(L)STTO(L*1) » SURF(L)) 
836 CONTINUE 
РО 839 (=1»+1М 
839 ADD2=ADD2+SURFIL) 
THIS SUB-PROGRAM DETERMINES THE TOTAL FLUX FOR OZONE, 
848 AREA2=AREA2 + ADD2 +AREA22 
WRITE (51»411) ADD2»AREA22,AREA2 
411 FORMAT(1X»>8H ADD2 = El7+.9»5X»10H AREA22 = El7+.995X»14HTOTAL AREA2 
l= €17.9) 
MUMP=MUMP-1 
ТӘ-ТЕМРЗ(КС) -273. 
GO ТО 830 
815 WRITE(51»842) TU 
842 FORMAT( 10H TF CO2 = F5e093X922H IS OUT OF TABLE RANGE) 


51 


883 


851 
854 
857 


412 


480 
450 


451 


410 
474 
490 
491 
413 


KN=KN-1 

TU=TEMP3(KN) -2736¢ 

GO TO 830 

CALL FORGEI(TUSTTCCIT)SRCCIT) ›, ВЕСО (КМ) »АКЕАЗ 3 } 

ADD3=0, 

IM=IT-1 

DO 851 L=leIM 

CALL FORGE2Z(RC(L) »eRC(L4+1) eTTC(L) oT TC(L4+1)eTURF(L) ) 

CONTINUE 

DO 854 L=lsIM 

ADD3=ADD3+TURFIL) 

THIS SUB-PROGRAM DETERMINES THE TOTAL FLUX FOR CARBON DIOXIDE. 
AREA3=AREA3+ADD3+AREA33 

WRITE (519412) АООЗ»АКЕАЗЗ»АКЕАЗ 

ЕОНМАТ( 1Х»8Н АОО2 - Е17.995Х%910Н АКЕАЗЗ - Е17.9»95Х» 14НТОТАІ АКЕАЗ 


l= Е17.9) 


THIS 598-РЕОСКАМ DETERMINES THE CO2 AND OZONE OVERLAP CORRECTIONS, 
JIP=3 

DO 450 1=1»20 

WU=CFW(KT) 

IF(WU-UC(1)) 450,470,480 

IF(WU-UC(I+1)) 470,473,450 

CONTINUE 

WRITE(51»451)WU 

FORMAT(17H LOG UF TAUCO2 - Ғ9,4»3Х»22Н IS OUT OF TABLE RANGE//) 
RCOLAPzxO, 

GO TO 452 

Iw=I 

GO TO (49024912492 )JIP 

CONTINUE 

CONTINUE 

Iw=I + 1 

60 TO 474 


4920CALL AOVLAP (WUSUCCIW)*sUCCIWEF1) $COXMISCIW)9*COXMISCIW*1)9COTAUF (KT) » 


452 


1ADD3»AREA33»DD»*RCOLAP ) 


JIP=JIP-1 
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DO 886 I=KOUNT ,+KT 
886 OG3TAUF(I) = EXPF(-0.1167 * SFW(1)) 

RO3LAP=O3TAUF (KT) * (ADD2+AREA22+GG) 

KIP=3 

КУ = 

DO 549 J=KOUNT | KT 

WU=CFW(J) 

DO 550 Iz1520 

IF (WU-UC(1I)) 55095702580 
580 IF(WU-UC(I*1)0) 57095739550 
550 CONTINUE 

WRITE(51+453) KASE(J) » WU 
4530FORMAT(1X»7H CASE =14,3X,16H LOG UCO2 TAU = F9e4,22H IS OUT OF TAB 

1ЕЕ RANGE) 

60 TO 549 
570 IW = I 
574 GO TO (59095915592) KIP 
590 CONTINUE 
573 IW=I+1 

GO TO 574 
5920CALL BOVLAP(WUSUC(IW) sUCCIW*1) s>COXMIS( IW) sCOXMIS( IW41) > 

1COTAUF(CJ) »5DUTCJ)sTCJ) 9TCJ*1) 5 CO2LAP( J) $DUT(J*1)0) 

KV=KV+1 

CO2LAP(KV) = CO2LAP(J) 
549 CONTINUE 

KIP = KIP - 1 

TCOLAP = RCOLAP 

ро 449 Ј=1, КҮ 
449 TCOLAP = TCOLAP + CO2LAP(J) 

DO 454 J=KOUNTsKT 
591 CALL BBVLAP(OSTAUF(CJ) GUTCJ) »GUTCJ. 1)9»3TCJ)sT(J*1) 903RLAP(J)) 
454 CONTINUE 

ТОЗКАР = RO3LAP 

DO 569 J=KOUNT KT 
569 TO3LAP = TO3LAP * O3RLAP(J) 

THIS SUB-PROGRAM DETERMINES THE TOTAL DOWNWARD FLUX FOR EVERY 
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C 


SOUNDINGe 
CALL ASUMF(LAREAL »AREA2 s AREA3sFF »6G»DD* TCOLAP +» ТОЗАР» ТОТА) 
WRITE(51979) (AREA1 sAREA2 »AREA3 9 FF 9GGsDD0s TCOLAP s TO3LAP sTOTAL) 


TIOFORMAT (2X sSHAREA] 9 3X 9 SHAREA2 9 3X 9 SHAREA3 9 6X s2HFF » 7 Х +» 22860» 6Х » 2800 » 


15X»6HTCOLAP»2X»6HTO3LAP»2X» 10HTOTAL FLUX//(1X»7E1749)) 


ГИУ = AREAl1+FF-TCOLAP/2+.-TO3LAP/2e 

FO3 = AREA2+GG-TO3LAP/2. 

ЕСО2 = AREA3+DD-TCOLAP/2e 

THE FOLLOWING VALUES HAVE UNITS CAL/CM SQUARE/DAY 

WRITE (519201) AREA] 9 FF sAREA2 »GGeAREA3 » DD oe TCOLAP s TO3LAP sFWVsFO3 >» 


1FCO2sTOTAL 


2010FORMAT(1X»10H AREAL = El7.9//1X1 H FF = E17.9//1X»10H AREA2 


1911 


1001 


- E17.9//1X5»10H GG = El7.9//1X»10H AREA3 - Е17.9//1Х,10Н DD 
= El7.9//1X»10H TCOLAP = E17e¢9//1X910H TO3LAP = El7.9//1X»10H 


3FWV - E17.9//1X»10H ҒОЗ = El7e¢9//1X910H FCO2 - Е17.9///1Х» 
423H TOTAL DOWNWARD FLUX = E17,9) 


XAREA1 = AREA1*CONST 

AFF = FF*CONST 

XAREA2 = AREA2*CONST 

XGG = GG*CONST 

XAREA3 = AREA3*CONST 

XDD * DD*CONST 

XCOLAP = TCOLAP*CONST 

XO3LAP - TO3LAP*CONST 

ХЕМУ = FWV*CONST 

ХЕОЗ = FO3*CONST 

XFCO2 = FCO2*CONST 

XTOTAL = TOTAL*CONST 

THE FOLLOWING VALUES HAVE UNITS WATTS/METERS SQUARE 
WRITE(C51»1011) KASE(KT) 

FORMAT(1Xs7HCASE = 14) 

wRITE(51»1001) XAREAL»XFF»XAREA2»XGG»XAREA3»XDD»XCOLAP»XO3LAP»XFWV 


1»XF03»XFCO2 »XTOTAL 


FORMAT(t 1X» 10H NAREAl = El7.9//1X»1 H NFF - Е17.9//1Х,10Н МАКЕА2 
- Е17.9//1Х,10Н МСС = E17,.9//1X»10H NAREA3 - E1749//1X»10H NDD 
= El7.9//1X»10H NCOLAP = El7+.9//1X»10H NO3LAP = El7+.9//1X>10H 
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229 


3NFWV =r eos 7 UN SOT! NOs 
423H TOTAL DOWNWARD FLUX = E17.9///) 


INTER=INTER-1 
KUT = KUT + 1 


МКАЗЕ (КИТ) = КАЗЕ(КТ) 
ZCOLAP(KUT) = XCOLAP 
ZO3LAP(KUT) = XO3LAP 


ZFWV(KUT) = XFWV 
ZFO3(KUT) = XFO3 
ZFCO2(KUT) = XFCO2 
ZTOTAL(KUT) = XTOTAL 
WFWV(KUT) с ХАКЕА1 + XFF 
WFO3(KUT) = XAREA2 + XGG 
WFCO2(KUT) = XAREA3 + XDD 


z El7.9//1X»10H NFCO2 


IF (INTER) 88859999777 
WRITE (519100) 

YKUT = KUT 

ADD1 = 0.0 

ADD2 = 0,0 

ADD3 = • О 

АМО11 = 060 

AVG22 - 0.0 

AVG33 - 0.0 

AVG44 = 0.0 

AVG55 = 0.0 

AVG66 = 0,0 

AVG77 = 0.0 

AVG88 = 0.0 

AVG99 = 0.0 

DO 8050 [=1 «KUT 

AVG11 = AVG11 + ZCOLAP(I) 
AVG22 = AVG22 + ZO3LAP (1) 
AVG33 = AVG33 -* ZFWV(I) 
AVG44 = AVG44 + WFWV(I) 
AVG55 = AVG55 + ZFO3(1) 
AVG66 = AVG66 + МРОЗ(1) 
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= Е17.9//1Х» 


AVG77 AVGI7 +°ZFCO2( 1) 

AVG38 AVG88 + WFCO2(I) 

AVG99 2 AVG99 * ZTOTAL(1I) 
8050 CONTINUE 

AVG] = AVG11/YKUT 


AVG2 = AVG22/YKUT 

AVG3 = AVG33/YKUT 

AVG4 = AVG44/YKUT 

AVG5 = AVGS5/YKUT 

AVG6 = AVG66/YKUT 

AVG7 = AVG77/YKUT 

AVG8 = AVG88/YKUT 

AVG9 = AVG99/YKUT 

DO 1111 I=1»KUT 

КАТТО1(1) = ZFWV(1)/ZTOTAL(I) 
RATIO2(1) = 2FO3(1)/ZTOTAL (1) 
RATIO3(1) = 22СО2(1)/2Т0ТА (1) 


1111 CONTINUE 
BO 1112 I=1»KUT 
ADD1 ADD1 + RATIOI(I) 
ADD2 ADD2 + RATIO2(1) 
ADD3 = ADD3 + RATIO3(1) 
1112 CONTINUE 


ADD11 = (ADD1/YKUT)*1002+ 
А0022 = (А0Ор2/ҮКЧТ ) +100. 
ADD33 = (ADO3/YKUT)*100, 


WRITE(5191022) (MKASE( 1) sZCOLAP(1) sZO3LAP(1) о МЕМУ ( 1) о ХЕМУ (1) » 
]ИРОЗ (1) »ДЕО03(1) »МЕСО2 (1) » 2 ЕСО2 (1) » ХТОТА (1) » 1=1»КкЏТ ) 


1022 ҒОНМАТ(1Х»85НСА5Е Н20-С02 Н20-03 UNCOR- СОККЕСТ- UNCOR- CO 
1ККЕСТ- | UNCOR- CORRECT- TOTAL /1X»88H OVERLAP OVERLAP REC 
2TED ED яу RECTED ED 03 RECTED ED CO2 DOWNWARD/1X»8 
34H WV FLUX FLUX 03 FLUX FLUX CO2 
4 FLUX FLUX FLUX/(1X91492F8039F100393F9e393F1003)) 


WRITE(5191031) AVG1sAVSG2 sAVG3 5 AVG4 sAVG5 sAVG6 9 AVG7 9AVG8 2 AVGI 
1031 FORMAT(IX»s4HMEAN/1X9»5HVALUESF7,3»F84.3»5F104, 3» 3F9, 3»3F10,4 3) 
WRITE(51,1022) ADD11,ADD22,ADD33 
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С 
С 
С 


1032 FORMAT(t1X»31HAVG VALUE CONTRIBUTION ON TOTAL/1X»32H DOWNWARD FLUX 
1 1М РЕКСЕМТ ве со осе »2-666 2 + 21862,» 2202) 
688 cND 


THIS SUBROUTINE EXECUTES FORWARD AND BACKWARD INTERPOLATION TO OBTAIN THE 
SPECIFIC ELASSER WATER R VALUE FOR EVERY REDUCED OPTICAL PATH/TEMPERA- 
ТОКЕ PAIR. 

SUBROUTINE FFIN(U9U1sU2sTsTleT2sFUIT1sFULT2sFU2T1 sFU2T2 oF) 
FFUTI=FULT1+(FULT2Z-FU1T1) *ABSF(CT-T1)/(T2-T1)) 
FFUT2sFU2TI1*(FU2T2-FUITI1)*ABSF(UOT-T1)/(CT2-T1)0) 
F1=FFUT1+(FFUT2-FFUT1)*ABSF ( (U-U1)/(U2-U1)) 

FFUT3 = FUIT2-(FUL1T2-FUITI)*ABSF((T2-T)/(T2-T1)) 
FFUT^-FU2T2-(FU2T2-FU2TI1)*ABSF((T2 T)/(CT2-T1)) 

F2= FFUT3+(FFUT4-FFUT3)*ABSF ( (U-U1)/(U2-U1)) 

Е-(Ғ1 ж» Ғ2)/2» 

END 


OSUBROUTINE FFOZ(USsU1s»U25» Ts T1» T2» GUIT1»*GUI1T2»GU2T1sGU2T2»G) 


GGUT1 = С6М1Т1%(60172-601Т1) ҒАВ5Ғ((Т-Т1)/(Т72-Т1)) 
GGUT2 - GU2T1«(GU2T2-GU2T1)*ABSF(CT- T1) /( T2- T1)) 
G1=GGUT1+(GGUT 2-GGUT 1) *ABSF ( (U-U1)/(U2-U1)) 
GGUT3=GU1T2-(GU1T2-GU1T1) *ABSF((T2 T)/(T2-T1)) 
GGUT4=GU2T2-(GU2T2-GU2T1) *ABSF((T2 T)/(T2-T1)) 
G2-2GGUT3* (GGUT4-GGUT3) *ABSF ( (U-U1)/(U2-U1)) 
G-(G1*G2)/2« 

END 


QSUBROUTINE FFCO2(U»U1sU2» Ts T1, T2»2DUITI1*3DUI1T25sDU2T1»DU2T2»0) 


DDUT1sDUITI*(DUIT2-DUITI)*ABSF((T-T1)/(T2-T1)) 
DDUT2 = DU2T1 + (DU2T2 - DU2T1) * ABSF((T-T1) / (T2 - T1)) 
D1=DDUT1+(DDUT2-DDUT1)*ABSF ( (U-U1)/(U2-U1)) 
DDUT3zDUIT2-(DUIT2-DUIT1)*ABSF((T2 T)/CT2-T1)) 
DDUT4= DU2T2-(CDU2T2-DU2TI1)*ABSF((T2-T)/(T2-T1)) 
D2 =DDUT3+(DDUT4-DDUT3)*ABSF ( (U-U1)/(U2-U1)) 
D=(D1+D2)/2e 

END 

SUBROUTINE RROZ(U9U1 9U2 sROU] sROU2 GG) 
GG=ROU1+(ROU2-ROU1)*ABSF ( (U-U1)/(U2-U1)) 

END 
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SUBROUTINE RRWAT(U»sULIsU2sRWU1 RWU2sFF) 
FFZRWU1*(U/URWU2-RWUI)*ABSFCCU-U1 )/(U2-U1) ) 
END 
SUBROUTINE RRCO2Z2(U 9U1] sU2sRCUL,RCU2 » DD) 
DDzRCUI*(RCU2-RCUI ) *ABSF( (U-U1)/(U2-U1) ) 
END 
SUBROUTINE ASUMF(AsB»CsD» EsF»X5, Y»sZ) 
Z = A+B+C+D+E+F-X-Y 
ЕМО 
ЗОВКОЈТТМЕ ВОУСАР(О 01» 02 оК] ои2»22 »01»11+12» ХХ»02) 
22= К] * (R2-R1)*A8GSF((U-U1)/(U2-U1)) 
XX=ZZ*(D1402)/2e*(T1I-T2) 
END 
SUBROUTINE AOVLAP(UsU1 9U292R19R29Z29A9B9DD 9X ) 
Z=R1+(R2—-R1)*#ABSF ( (U-U1)/(U2-U1)) 
X=Z*(A+84DD) 
END 
SUBROUTINE BBVLAP(AsR1 s9R20T1l9T29Z) 
Z= A® ((R1 + R2)/2¢ * (T1-T2)) 
END 
SUBROUTINE FFFN(UsU1 9U25R1 sR25XX) 
AX=R1+(R2—-R1) *ABSF((U-U1)/(U2-U1)) 
END 
SUBROUTINE FORGEI1(TFsT1lsR1sRFyx) 
X=(R1+RF)/2e*%(TF-T1) 
END 
SUBROUTINE FORGE2(RI»R2»T1»T2,Y) 
ҮсАВОҒ((К1%К2)/2.%(Т2-Т1)) 
END 
END 
FINIS 
-ЕХЕСОТЕ» 
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Appendix I B. Modification Instructions to Program RADFLUX (FORTRAN 63) 


1. Each sounding is subdivided into levels starting from the interface 
at level I = 1 to the top of the atmosphere (0.1 mb) at level I = KZ. 
KZ is the number of the last level in the sounding. 

2. The following must be completed in order to determine F, from any 
preselected reference level different from the interface level: 

a. Note the count or number of the pre-selected level in the sound- 
ing and introduce this number as shown below in the example. 

1. Suppose selected reference level is 7 out of 30 levels in 
a sounding, hence, punch this number on a data card with 
the unit digit in column 4. This reference level card is 
the second data card in each of the case atmospheres. Do 
this procedure anytime the reference level changes for any 
case-atmosphere. 
3. If a level is added or subtracted in a sounding listing, the number 
on the first data card of the respective case atmosphere must be changed 
accordingly. For example: 

a. Suppose that there are 29 levels in a sounding, thus, the first 
data card in the sounding must read 29; however, suppose that 1l 
level is added to the sounding, thence, the first data card must 
be changed to read 30. 

b. One must place the unit digit of the specific number in column 
4 of the data card. 

4. Finally, the INTER (constant) card must replaced when the number 
of case atmospheres in the study changes. For example: 
a. Presently there are 63 case atmospheres as data input to program 


RADFLUX, thus, the INTER card reads "INTER = 63." 
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Ы, Suppose the data input changes to 100 case atmospheres; thence, 
the INTER card must be changed to read "INTER = 100." 
5. Program RADFLUX was written for FORTRAN 63, however, this program 
can be made operational in any FORTRAN computer language provided 
that minor modifications to the control cards are made. One must 
check with the respective computer language manuals for further 


information. 
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Appendix I, Section C. Computer Program SUMFLUX (FORTRAN 63) 


-COOP,BOX T» TUPAZ J»S/1S/28»20»20000. 
-FTN»*L»Ee 


C 
С 
e 


PROGRAM SUMFLUX 
THE PURPOSE OF THIS PROGRAM IS TO INTEGRATE RADIATION INTENSITY 


CASE» 
BETWEEN O AND 90 DEGREES TO OBTAIN TOTAL UPWARD FLUX FOR EVERY 
ODIMENSION KASE(7) »PRES(7)» THETAU7) X07) X2(73 3UH2(7) 3UO3(7) 5 X40 7)» 
LTEMP (7) sSTOTINT(7) sCHGSINT(7) 
DIMENSION KOUNT( 100) »sPRESS(100) 3 TEMPER( 100) 5 CHFLUX( 100) » 
1BBFLUX (100) 
INTER = 62 
KUT = 0 
CONST=0.56687E-7 
PI = 3.14159 i 
RADCON = 0.017453 
KZ = 6 
ТНЕТА(6) = 90,0 
7770READ (50»2) KASE(1)»PRES(1)»>THETA(1)>X(1)>»X2(1)»UH2(1)»U03(1)»Xx4(1 
1}, ТЕМР (1) 
2 FORMAT(13»F9.1>»F8.0>»F10.3»2F7.3»F6+.3»F10+.3»F701) 
KT = KZ-1 
OREAD (5053) (KASE(I) sTHETACI) »X( 1) sX2¢1) oX4(1)91=29KT) 
3 ЕОКНМАТ(13»9Х», 2860» 210, 3,27 63»13Х» Е1063) ; 
SUM1=00 
SUM2=0 • 
DO 100 Is1»KT 
CALL FORD1UX(3)2)»XCA) 93X(5) »X6( 3) »X&(4) »X64(5)»X(6) »X6(6)) 
CALL FORD2(X( 1) oXC 141) oX4(01) X40 141) sTHETACI) + ТНЕТА (1+1) + ТОТІМТ (І) 
1+CH4INT(I)»RADCON) 
WRITE(51»5) KASE(I)»I»TOTINT(])>»CH4INT(I) 
SUFORMAT(1X»8H CASE = I5»5X»5H ] = I4»5X»10H TOTINT = El7.9>»5X»10H С 
LHGINT = E17e9//) 
SUM1=SUM1+TOTINT(1) 
SUM2=SUM2+CH4 INT(I) 
100 CONTINUE 
WRITE(5194) (КАБЗЕ(1)»РВЕЗ(1)»ТЕМР(1) + (1»ТНЕТА(1 )» Х(1)»Х4 (1)»1=]»К2 
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1)) 
40FORMAT(1X»8H CASE * I55»5X»12H PRESSURE = F5.0,5X,8H TEMP = F5.0// 
1(1Xs5H І = I4»5X»10H THETA1 = F5.1»5X»,19H TOTAL INTENSITY * F6,23» 
25X»18H CH 4 INTENSITY = Рбез//)) 
FLUX]=SUMIXP ] 
FLUX2=SUM2#P] 
wRITE(51»6) KASECI)»FLUXISFLUX2 
6 FORMAT(1X»8H CASE * I555X»9H FLUX1] = E17,9»5X5»9H FLUX2 с Е17.9///) 
CALL FORD3( CONST» TEMP(1) »FLUX3) 
WRITE(51»300) FLUX3 
300 FORMAT(1X»9H FLUX3 с Е17.9///) 
PRESLN = LOGF(PRES(1)) 
INTER = INTER-1 
KUT = KUT + 1 
KOUNT( KUT) = KASE(1) 
PRESS(KUT) = PRES(1) 
TEMPER( KUT) = TEMP(1) 
CHFLUX(KUT) = FLUX2 
BBFLUX(KUT) = FLUX3 
IF (INTER) 88899999777 
999 wRITE(51»200) 
200 FORMAT (1H1) 
WRITE(519400) (KOUNT( I) sPRESS(])»sTEMPER( I) »sCHFLUX(!) ЭВВЕІ ӘХ(І)» ігі 


l»KUT) 

400 FORMAT(1X»53HCASE INTERFACE INTERFACE CHANNEL 4 BLACK 
1/1Х»52Ң PRESSURE TEMPERATURE FLUX BODY/49X» 
2ӘНҒЦОХ//(1Х»14»Е8.0»Ғ13.0»Ғ15,3»Ғ16.3)) 

888 ЕМО 


SUBROUTINE FORDI1(As»sB»sC»D»EsF »Gs H) 

С = С(+0.6866*%А -1.8649*B +1.1783*C 

Н = Е+0,6866%0 -1.8649*E +1.1783*F 

END 

SUBROUTINE ҒОКО2(А»В»С»ОЭЕ»Ғ»Х,Ү92) 
X=(A+B)/2+* ((SINF(IF*Z))**2 - (SINF(EXZ))**2) 

Y = (C+D)/20*((SINFIFAZ))*%2 — (SINF(EXZ))**2) 
END 
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